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O ! Abstract 



We study the collision of two Reissner-Nordstrom gravitational shock waves in AdS and 
show that the charge completely prevents the formation of marginally trapped surfaces of the 
Penrose type with topology S^~^, independently of the energy and the value of the impact 
parameter. In the case of head-on collisions, a numerical analysis shows that no trapped 



H ' surfaces with topology S^ x S^ ^ form either. 
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Introduction. The study of collision of gravitational waves has been a very active area 
of research in the field of general relativity (see pLj for a summary of results). The interest 
in the subject has been renewed in the context of the AdS/CFT correspondence [2j. The 
underlying reason is the expectation that colliding waves in AdS spacetime could provide a 
reliable gravitational dual of the high-energy collision of two "nuclei" , modeled by energy lumps 
in the holographic strongly coupled gauge theory [3l H]. Although these energy lumps do not 
reproduce all the properties of heavy ions, these theoretical experiments can be used as models 
where collective gauge theory effects, relevant for realistic heavy ion collisions, are tractable. 
An expected outcome of the collisions of two gravitational waves is the formation of a black 
hole. Holographically, black hole formation can be interpreted as the thermalization of the 
gauge theory plasma resulting from the collision of the lumps. 

In spite of the recent progress in the numerical study of the problem of the collision of two 
gravitational waves [5], a full understanding of the highly nonlinear dynamics dominating the 
physics in the interaction region is still missing. An alternative to a full numerical simulation 
of the collision process (in the spirit of numerical general relativity) is to look for the formation 
of a trapped surface that would signal the eventual presence of an event horizon [3] . Although 
in many cases this also involves numerical analysis, the approach is technically simpler. The 
idea [6j consists in looking for marginally trapped surfaces with topology S^^"^, lying on the 
hypersurface {u = 0,v < 0}U{u < 0,v = 0}, which does not involve solving for the geometry in 
the interaction region. The information about the collision of the waves is nevertheless encoded 
in the nontrivial matching condition aX u = v = 0. In the following, we will refer to these 
marginally trapped surfaces as of Penrose type. 

Varying the parameters of the collision, it is possible to find thresholds for the formation 
of Penrose trapped surfaces. Holographically, this can be seen as a threshold for the onset of 
the plasma thermalization after the collision takes place. This is what happens, for example, 
for large enough impact parameter [3, E], or when the sizes of the two colliding lumps are very 
different [8]. A more intriguing threshold was found in [H] associated with the transverse size 
of the gravitational source of the shock wave in AdS, which however does not induce a change 
in the holographic energy-momentum tensor. 

Trapped surfaces with topology 5^ '^ in the collision of RN-AdS shock waves. It 

would be desirable to extend the analysis of the formation of marginally trapped surfaces of 



Penrose type to more general types of incoming waves, in particular those obtained taking an 
infinite boost limit [TU] of a D-dimensional Reissner-Nordstrom (RN) solution asymptotically 
anti-de Sitter (AdS) 
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where L is the radius of AdS, and M and Q are given in terms of the mass m and electric 
charge q by 
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with fijv the volume of the A^-dimensional sphere S^ . Performing a boost with Lorentz param- 
eter 7, the shock wave geometry is obtained in the limit 7 — ?■ cxo while keeping 



fi = jM, 
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fixed [m [T2I [13]. The metric takes the form (in Poincare coordinates) 
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The metric function $(q')RN depends on the AdS chordal coordinate 
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and can be expressed as 



$(g)RN = ^q)D - ^$(g)2D-3, 
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where ^{q)D is the profile of the Aichelburg-Sexl shock wave in AdS^i [151 E]- Therefore we 
have 
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To avoid confusion, we recall that the energy parameter fi = 7M used here and in [T3] is related 
to the parameter E = ■ym of Ref. ^ by the rescaling E = ^^^Q£)_2fJ'- 

We analyze the problem of collision of two shock waves. Generically, in the region outside 
the interaction wedge, the metric reads 
L2 
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It is convenient to change coordinates (m, f , z, xt) =^ {U, V, Z, Xt) to remove the distributional 
terms in the metric of the incoming waves (see the Appendix of Ref. PJ). Following Penrose 
[6], we look for surfaces S with topology S^~'^ contained in the past ligh-cone {U < 0, V^ = 
0} U {U = 0, V < 0}. This is parametrized in terms of two nonnegative functions ip±{Z, Xt) as 

U = -^+{Z,Xt) ^ = I ^ = ° 

V = Q ' ~ \ V = -^-{Z,Xt) 

where 5 = 5+ U iS_. Defining 
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^±(Z,Xt) = -V'±(^,^t) 
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the condition that the expansion of the congruence of null geodesies normal to S vanishes reads 
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where $-1- are the profiles of the two incoming waves. 

In the case of head-on collisions, we can exploit the 0{D — 2) invariance of the system and 
assume that '^±{Z,Xt) = ^±('?) only depends on the chordal coordinate (^. In terms of it, 
the metric of the hyperbolic transverse space takes the form 






and the trapped surface equation ( TTTT) is given by 
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Since S has two branches, we have to make sure that the null geodesies are continuous 
across the (D — 3)-dimensional surface C = Sj^r\S^ defined by \E'±(g) = 0. This is implemented 
by the boundary condition 



5«U*±96*= 



(14) 



where Qab is the metric on H[£)_2. Using the chordal coordinate, the surface C is parametrized 
hj q = qc, and the previous condition can be written as 

4r2 or 

The sign is fixed by the requirement that the trapped surface has topology S^^"^, which means 
that \E'±(q') > for < q < qc- The trapped surface is determined by the solution to the 
differential equation ( TT3l) subjected to the two boundary conditions 

^±(gc) = 

^'±fe) = - , f^ , ^, (16) 

^/qc{qc + ^) 

For simplicity, let us assume collision between two identical shock waves, so we have '^+{q) = 
^_(g) = "^{q) and $+(g) = $-(q') = '^(q')RN- The general solution to the differential equation 
( TT3l) can be written 

^(g) = '^(g)RN + C7i$i(g) + C2$2(g), (17) 

where Ci^2 are two integration constants and $1,2(0') are the two independent solutions to the 
homogeneous differential equation 

$i(g) = l + 2q 

Uq) = q'-''2F^(^D-2,^;D;-^y (18) 

If the trapped surface has topology S^^'^, the function \l/(g) has to be regular for < g < gc 
and this forces C2 = 0. Imposing besides $i(gc) = determines the remaining constant to be 

Ci = -— ^$(gc)RN. (19) 

l + 2qc 
The value of qc is now determined by the second condition (fT6|) . Plugging 

*(g) = $(g)RN - -^^|^$(gc)RN, (20) 



into the second equation in (]T6l) . we find the algebraic equation 
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Figure 1: Plot of the function \l/(g) for D = A (left panel) and D = b (right panel). In all cases 
the energy of the head-on collision is Gat/x/L^"^ = 1 and the charge parameter (from top to 
bottom) VG^e/L^-^ = 0.5, 0.75 and 1.0. 



The problem of the existence of trapped surfaces of the Penrose type in the collision of two 
RN-AdS shock waves has been studied in [13] (see also [13] ), as well as in a first version of 
this paper. The analysis presented there, however, has a fundamental flaw consisting in that 
the function ^{q) has a second zero below qc and therefore the surface defined by it does not 
have topology 5*^"^ as assumeqj. The origin of this problem lies in the fact that the charge- 
dependent term in $rn(q') shown in Eq. ([7]) tends to minus infinity as g — )■ 0^. Indeed, whereas 

4— D 

in this limit the /i-dependent term diverges as q~^^ for D > 4 and — log q for D = 4, the e- 

7— 2D 

dependent term diverges as —q~^ for D >A. This second term dominates near g = for any 
e > (and i5 > 4), so \E'(g) goes to minus infinity as g — )■ O"*" [see Eq. (120|) ]. This behavior is 
explicitly shown in Fig. [H where the data for the function \E'(g) is plotted for D = A and D = 5 
and various values of the charge. 

On physical grounds it is to expect that this problem remains when considering a non- 
vanishing impact parameter for the collision. To check this requires solving numerically the 
Laplace- type equation f lTT]) with the appropriate boundary conditions [161 El E] • We work in 
radial coordinates r = 2Lyg(g + l), where the metric of the transverse hyperbolic space takes 
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The system is 0{D — 3)-invariant and the Penrose trapped surface is parametrized by \l/±(r, 9). 
To solve for these functions, it is convenient to define Hj-{r, 9) = $±(r, 9)^^ — \l/±(r, 9), where 
$-i-(r, 6')rn are the profiles of the incoming RN-AdS shock waves whose sources are located 



respectively at r-t = |, 6^+ = 0, 9- 
H±{r,9) satisfies the equation 
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This has to be solved in a domain < 9 < rr, < r{9) < LG{9). As a consequence of the 
symmetry of the problem, the functions $±(r, 6')rn and '$±{r,9) satisfy 

^±(r, 9) = vl/^(r, n-9), $±(r, 9)^^ = <l>^(r, vr - 9)^^. (24) 

This implies Neumann boundary conditions for the function Hj-{r,9) on the lower boundary, 
deH±{r,0)=0 = deH±{r,n). 

The function G{9) determining the shape of the candidate trapped surface also has the 
symmetry G{9) = G{'k — 9). It has to be chosen such that the following conditions are satisfied 
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r=LG{e) 

We solve this boundary problem for the collision of two identical shocks using the method 
devised in [16]: we solve Eq. ( l23l) numerically in a 50 (angular) x 100 (radial) grid and find G{9) 
through a trial-and-correction loop. In the following we summarize our results for the collision 
of two waves of the same energy. The details of the implementation of the method can be found 
in [8] , where it was applied to the off-center collision of two Aichelburg-Sexl-AdS shock waves 
in various dimensions. 

The numerical data obtained shows that \l'(r, ^) not only vanishes at the boundary r = 
LG{9), but also in the interior. In Fig. [2] we have plotted the radial profile of the function 
\E'(r, ^o) at ^0 = 0, both for D = 4 (left panel) and D = 5 (right panel). The conclusion is that 
there are no trapped surfaces of the Penrose type with topology S^^"^ produced in the collision 
of two RN-AdS shock waves, both with or without impact parameter. 
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Figure 2: Plot of the section ^ = of the function \I/(g, 6) in four (left panel) and five dimensions 
(right panel) with nonvanishing impact parameter b/L = 0.3. Again, the energy of the collision 
is Gi\ifi/L^^^ = 1 and the charge parameter (from top to bottom) a/G/vc/L^"^ = 0.5, 0.75 
and 1.0. 



Changing the topology. Given the results above, it seems natural to look for marginally 
trapped surfaces with topology S^ x S^~^. Mathematically, this means that we allow for the 
possibility of \l/(g) vanishing at two points that we denote by gin and gout, and such that 



^±(g) > for gin < g < gout- 



(26) 



The surface C has therefore two components Cin and Cout defined respectively by g = gin and 
Q = gout- Let us analyze now the case of a symmetric head-on collision. The condition f lT^ 
for the continuity of the congruence of normal null geodesies across C now translates into the 
couple of equations 
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The choice of signs is fixed by Eq. ( 126|) . Hence, we have to solve the differential equation ( 1T3|) 
with the boundary conditions 

(28) 

(29) 

2L , , 
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In the case of the solutions found in [13j, it is important to notice that Eq. fl3T]) is not 
satisfied at the internal zero of "^{q)- Hence, to have a chance of finding trapped surfaces 
of nontrivial topology we have to be more general. An important change introduced by the 
S^ X S^~^ topology of the trapped surface occurs in the form of the general solution to Eq. 
f lT5]) . Now the region of interest excludes both g = and g = oo, so there is no reason to set 
C2 = in (TT7|) as we did when assuming that the trapped surface had topology S^^"^. In fact, 
Ci and C2 are determined by Eqs. fl28|) and fl29|) 



(^l$l (gout) + <^2$2 (gout) = -$(gout)RN, 

Cl$l(gin)+C2$2(gin) = -$(gin)RN. (32) 

Once these constants are solved in terms of gin and gout) we impose the conditions (1301) and (13 ip . 
This provides two algebraic equations that, in principle, are enough to determine the values of 
the internal and external radii of the trapped surface. 

It is interesting to notice that the resulting function \E'(g) is independent of the mass pa- 
rameter fi. The right-hand side of both equations in fl32|) is the sum of two terms, respectively 
of order fx^ and /x. We write Ci^2 = C*i 2 + A^W 2; whereas Eq. ([6]) can be recast as 

$(g)RN = AD/i$2(g) + $(g) with An = L ^'' ^^^^^ {S^) ■ ^^^^ 

Solving the system ( l32i) order by order in /i gives C| = 0, C2 = —An- From Eq. ( TT71) we 
find that the //-dependent part in \E'(g) cancels. 

A numerical analysis of the system of equations fl30|) -f l3T]) . with the values of the constants 
found from ([32]), renders no nonvanishing solutions for g^n and gout- This implies that there are 
no trapped surface of topology S^ x S^~^ formed as the result of the head-on collision of two 
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Figure 3: Plot of the implicit equations (130|) and ( 13T]) for two values of the charge in D = 5. The 
two curves approach each other close to the origin but do not cross. The dashed line represents 
the diagonal gin = gout, and shows that both curves lie in the "physical" region gout > gin- 

RN-AdS shock waves. This is pictorially illustrated in Fig. |3l where the two curves defined by 
Eqs. fISUl) and fl^ do not cross each other outside the origin. 



Closing remarks The study of the problem of collisions of two RN-AdS shock waves pre- 
sented here shows that the charge parameter e completely prevents the formation of marginally 
trapped surfaces in the region {u = 0,f < 0} U {m < 0,f = 0} with topology S^~^, indepen- 
dently of the value of the impact parameter. In the case of head-on collisions, we have not 
found any trapped surface with topology S^ x S^~^ in the same region either. 

These results do not preclude the existence of trapped surfaces in other regions of the 
spacetime formed as the result of the collision. In [T71 [12] the problem of the collision of both 
Aichelburg-Sexl and RN shock waves in flat space was analyzed, looking for trapped surfaces 
lying in the future light-cone {u = 0,v > 0} U {u > 0,v = 0}. In the case of AdS, the analysis 
needed to find these trapped surfaces is more involved than the one required to find the Penrose 
trapped surface, since the change of coordinates needed to eliminate the distributional terms 
in the metric is nontrivial in the "future ligth-cone" region. As a consequence, the equations 
determining the trapped surface are much more complicated and have to be solved numerically, 
even in the case of head-on collisions. This problem is under current investigation. 
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